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Abstract 

The status is reviewed of recent high precision measurements of inclusive-jet, dijet and trijet production in deep- 
inelastic scattering and photoproduction by the HERA experiments HI and ZEUS. The measurements are in good 
agreement with perturbative QCD calculations at next-to-leading order and are used for the extraction of the value of 
the strong coupling at the mass of the Z boson, a s (Mz). The methods and results of the QCD analyses are presented 
and a summary of the values of a s (Mz) from recent jet measurements at HERA is given. 
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1. Introduction 

Since the advent of particle colliders the study of 
jet production has been essential for our understand- 
ing of the strong force and the development of quantum 
chromodynamics (QCD). Nowadays perturbative QCD 
(pQCD) is firmly established and jets are used exten- 
sively at the Tevatron and the Large Hadron Collider 
(LHC) in analyses searching for new physical phenom- 
ena. In these analyses an overwhelming amount of 
background originates from Standard Model (SM) QCD 
processes. Therefore potential discoveries depend heav- 
ily on our understanding of QCD, which can be im- 
proved by precision jet measurements. 

One of the key predictions of QCD is the running of 
the strong coupling a s as function of the renormalisa- 
tion scale /u r . However, the absolute value of a s at some 
starting scale has to be obtained from experimental data. 
The direct sensitivity of the jet production cross sections 
to a s makes jets an ideal tool for precision measure- 
ments of a s . Additionally, since jet data usually span 
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over a large range of values of fx r , it can be used for 
rigourous tests of the evolution of a s . 

In contrast to hadron-hadron collisions, in deep- 
inelastic scattering (DIS) there is only one proton 
present in the initial state. Therefore DIS provides a 
very clean environment for jet production. It is possi- 
ble to study pQCD without the complication of having 
to disentangle contributions from pile-up and multiple 
interactions, and only a minimum amount of modelling 
of non-perturbative effects is needed. 

Inclusive neutral current (NC) and charged current 
(CC) DIS data are used in analyses of the proton struc- 
ture for determinations of the parton distribution func- 
tions (PDFs) of the proton. These are one of the key 
ingredients for predictions at the LHC and uncertainties 
on the PDFs are reflected in the theoretical predictions. 
Inclusive DIS data are sensitive to the gluon density and 
a s only in next-to-leading order (NLO) through scaling 
violations, however. Jet data provide direct sensitivity to 
the gluon density and a s already in leading order (LO). 
The inclusion of DIS jet data reduces the correlation be- 
tween the gluon PDF and a s , thereby making a simul- 
taneous determination of proton PDFs and a s with DIS 
data alone possible 1 1 1. 

Jet production in photoproduction, where a proton 
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collides with a quasi real photon, additionally provides 
information on the structure of the photon. Since in 
these interactions the photon can show hadronic struc- 
ture, this kinematic regime closes the gap between DIS 
and hadron-hadron collisions. Having obtained the pro- 
ton PDFs using DIS data, jets in photoproduction can 
also be used for independent tests of the obtained pro- 
ton PDFs. Additionally, in photoproduction the specta- 
tor partons in the proton not participating in the primary 
interaction, may interact with partons of the resolved 
photon. In these processes multiple interaction models 
can be studied in an environment free of pile-up. 

The HERA collider provides a unique opportunity to 
study jet production in electron-proton (ep) collisions at 
the highest available centre-of-mass energies. After its 
shutdown in the year 2007 the HI and ZEUS collabo- 
rations invested a lot of effort in improving even further 
the detector calibrations and reconstruction algorithms. 
The benefits of these efforts become visible now with 
jet measurements of unprecedented precision. In this 
document the most recent precision jet measurements 
at HERA are reviewed, together with the QCD analyses 
and determinations of a s from jet cross sections. 

2. QCD and the running coupling 

One of the central predictions of QCD is the de- 
pendence of the strong coupling on the renormal- 
isation scale p r , which can be obtained from the 
renormalisation-group equations. It is expressed in 
terms of the /^-function 

P{a s (jUr)) = fit —5 , ( 1 ) 

which can be expanded in powers of a s and the result- 
ing coefficients /3„ are known up to four loops \2\. The 
behaviour of a s as function of the scale yU r is one of the 
best-known quantities in QCD. However, the absolute 
value of a s at some starting scale po has to be obtained 
from experiment. Usually the choice fio = M z is used, 
where M z is the mass of the Z-boson. Uncertainties 
in the value of a s (M z ) are reflected in uncertainties of 
SM predictions at the LHC such as W, Z, ft or Higgs 
production. The spread of the predicted cross sections, 
usually combined with the uncertainty from proton PDF 
parametrisations, can be as large as 10% [3|, showing 
the need for precision measurements of a s (M z ). Ad- 
ditionally, since one of the key predictions of QCD is 
the universality of a s {M z ), testing the compatibility of 
the obtained values of a s (M z ) from a variety of differ- 
ent processes measured at different scales p. t , provides a 
stringent test of QCD itself. 




Figure 1: Some important Feynman diagrams of neutral current ep- 
scattering. The LO process is shown on the left, followed by one in 
each case of the boson-gluon fusion and QCD Compton diagrams. 

3. Deep-Inelastic Scattering 

The LO Feynman diagram of ep scattering is shown 
in Figure [T] (left). The incoming electron with four- 
momentum k interacts with the proton with four- 
momentum P through the exchange of a virtual boson 
with four-momentum q — k-k' . In NC reactions, where 
the exchanged boson is either a photon or a Z-boson, 
the final state consists of a scattered electron with four- 
momentum k' and the hadronic final state X. 

The kinematics of the scattering process can be 
described by the negative four-momentum transfer 
squared Q 2 - -q 2 , also called the virtuality of the 
exchanged photon, the Bjorken-scaling variable x = 
Q 2 /(2P ■ q) and the inelasticity y = (P ■ q)/(P ■ k). These 
variables are related at fixed centre-of-mass energy -^s, 
through the relation Q 2 = sxy, such that two variables 
are sufficient to uniquely define the scattering process. 

The variable Q 2 is also used to separate the ep scat- 
tering process into two different kinematic regimes. At 
vanishing Q 2 the proton collides with a quasi-real pho- 
ton and this y*p interaction is referred to as photopro- 
duction. For Q 2 significantly larger than 1 GeV 2 and the 
invariant mass of X much larger than the proton mass, 
one speaks of DIS. 

The inclusive NC scattering cross section for the re- 
action e ± + p — > e ± + X can be calculated in terms of 
the structure functions F2, xFj, and Fi (4). The struc- 
ture function F2 gives the dominating contribution to 
the cross section in the bulk of the phase space. It 
is in LO related to the quark densities of the proton, 
F2 = xYuie 2 {qi + qi), where the sum runs over all ac- 
tive quark flavours, and the charge <?,■ of quark i is given 
in units of the elementary charge e. The structure func- 
tion xFj is directly proportional to the difference be- 
tween the quark and anti-quark densities, and Fl is sen- 
sitive to the gluon density. However, contributions to 
the cross section from the gluon density, either directly 
through Fl or because of corrections to F2 (scaling vi- 
olations), are non-zero only at NLO and beyond. Fur- 
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thermore, at NLO the gluon density enters the DGLAP 
equations multiplied by a s , which leads to a large cor- 
relation between the gluon density and a s in determina- 
tions of PDFs of the proton. 

4. Jet production at HERA 

While the inclusive ep scattering process is at LO 
purely of electroweak nature, and QCD enters only 
through higher-order corrections, in jet production QCD 
plays a dominant role already at LO. This leads to a 
direct sensitivity of the jet cross sections to the strong 
coupling a s and the gluon density of the proton. Hence, 
the study of jet production provides a prominent way 
to determine the value of a s (Mz) as well as important 
complementary information for QCD analyses of struc- 
ture function data . 

Jets are collimated sprays of hadrons, which are 
created by the showering and hadronisation of outgo- 
ing final state partons. The exact definition of a jet 
depends on the employed jet algorithm, which speci- 
fies which hadrons are clustered into a jet and how the 
jet four-momentum is calculated. At HERA, usually 
the longitudinally invariant kj algorithm J6] |7] with the 
P T recombination scheme |8| is used. Recently also 
studies with the anti-£ T and the SISCone [ 10 1 algo- 
rithms have been performed ifTTl . In DIS analyses the 
jets are clustered in the Breit frame of reference, where 
the transverse momentum Pj of jets is a measure of the 
hardness of the underlying QCD process. In analyses of 
jet production in y*p interactions, the jets are clustered 
in the laboratory rest frame. 

4.1. Jet production in DIS 

The cross section for jet production in the Breit frame 
in DIS can be written as a series expansion in powers of 

o> = J] \<%filp(£>ti) (2) 

where ff/ p are the PDFs of the proton, and the coeffi- 

(ft) 

cient functions C. can be calculated in pQCD up to 
some order of n. The variables pf and p r are the fac- 
torisation and renormalisation scales. The last term in 
Eq. (|2} represents the hadronisation correction, which 
is due to non-perturbative effects from the hadronisation 
process. The two Feynman diagrams in the centre and 
on the right hand side of Fig. [T]show two of the LO, i.e. 
0(a s ), contributions to the jet cross section. The boson 



gluon fusion (BGF) processes introduce a direct sensi- 
tivity of the cross section to the gluon density of the 
proton and dominate the cross section at small values 
of Q 2 and transverse jet momentum P T . With increas- 
ing values of Q 2 and P T , the QCD Compton (QCDC) 
processes become dominant. These are sensitive to a s 
and the quark distributions and can therefore help to re- 
duce the correlation between the gluon and a s in QCD 
analyses. The coefficient functions C w are known up to 
NLO for dijet and trijet production, resulting in 0(a 2 ) 
calculations for inclusive jet and dijet cross sections and 
0(al) calculations for trijets. The dominating theoreti- 
cal uncertainty is typically due to missing higher orders 
in the perturbative series. It is estimated by arbitrary 
variations of p r and pf by conventional factors of 0.5 
and 2, and is typically of the order of 10-15%. 

4.2. Jets in photoproduction 

In the case of photoproduction the photon can either 
interact directly with a parton in the proton or it can 
act as a source of partons, one of which interacts with 
a parton in the proton. In the latter case Eq. (|2]) has 
to be modified to account for the resolved photon con- 
tribution JT2). This is achieved by convoluting the co- 
efficient functions with the photon PDFs in addition to 
the convolution with the proton PDFs. The coefficient 
functions then become cffj (x, x y ,p 2 ,p 2 ^, where the in- 
dices i and j account for the interacting constituents of 
the proton and the photon. The variable x y is the longi- 
tudinal momentum fraction of the photon taken by the 
parton participating in the hard interaction. In LO dijet 
production, x y can be expressed as 

„ bs _ PT,ie-"' +PT, 2 e-»? 

where Pjj and 77, are the transverse momenta and pseu- 
dorapidities of the two hard jets, and E 7 is the energy of 
the exchanged photon. In the LO picture x° hs = 1 cor- 
responds to the point-like component of the scattering 
process. At NLO this identification no longer holds, and 
a sizeable contribution of the hadronic component can 
be observed IFPJI . In order to relate the j*p cross sec- 
tion to the total ep cross section, the photon flux is usu- 
ally calculated using the Weizsacker- Williams approxi- 
mation. 

The transition from photoproduction to DIS is inter- 
esting for a number of reasons. In jet production in y*p 
interactions the only hard scale is the jet Pj, whereas 
in DIS Q may become as hard or even harder than Pj. 
The presence of a second hard scale poses an interest- 
ing problem for the choice of the renormalisation scale 
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in pQCD calculations in DIS, which is of no concern 
in photoproduction. Furthermore, in DIS no underly- 
ing event is present, and jet production can be studied in 
a clean environment without contaminations from non- 
perturbative contributions. In photoproduction the con- 
stituents of the resolved photon which did not undergo a 
hard interaction may scatter from the spectator partons 
of the proton, such that multiple-interaction models can 
be studied, and information about the transition from 
lepton-hadron to hadron-hadron scattering may be ob- 
tained. Lastly, in a scenario that is analogous to the DIS 
case, the gluon density of the photon may be constrained 
by the study of jet production in y*p interactions. 

5. Determination of a s from jet cross sections 

The HI and ZEUS collaborations have adopted dif- 
ferent methods of extracting the strong coupling a s (Mz) 
from jet cross sections. In the following, both methods 
are briefly described and their advantages and disadvan- 
tages are summarised. 

5.1. Determination of a s (M z ) by HI 

The HI collaboration employs the Hessian 

method |[T4l IT31 for deriving values of a s (Mz), 
where the function 



X 



V T M- 1 



(4) 



is minimised. In this definition the matrix M is the 
full correlation matrix, which takes into account the 
uncorrelated statistical and systematic uncertainties as 
well as the statistical correlations between different data 
sets lfT6l . The vectors V describe the difference between 
the experimental data and the theoretical predictions 



T J eX P 1 

Vi = cr j - cr, 



,_theo 



1 - ^ A ik e k 



(5) 



where <x™ p and <x' heo are the experimental and predicted 
cross sections for a given bin i, and A,-* is the effect of 
the correlated systematic uncertainty k on the measured 
cross section in bin i. The Hessian parameters e k are 
left free in the fit and used as penalty terms. With this 
method the experimental uncertainty is given by the nat- 
ural choice of y 1 — v 2 . + 1 . 

A A mm 

The theoretical uncertainties due to missing higher 
orders in the NLO calculations, hadronisation correc- 
tions and the PDF uncertainties are determined with the 
offset method. In this procedure the theoretical predic- 
tions are recalculated with varied input parameters ac- 
cording to their uncertainty and the fit is repeated. The 



obtained relative change in a s (Mz) is stated as uncer- 
tainty. 

5.2. Determination ofa s (Mz) by ZEUS 

The ZEUS collaboration typically performs QCD 
analyses using NLO calculations with PDF sets that 
were obtained for different values of a s (Mz). The de- 
pendence of the theory cross section on a s (Mz) as func- 
tion of an observable X is parametrised in each measure- 
ment bin i by a second-order polynomial in a s (Mz), 



dcr 
dX 



c M a s (M z ) + c 2 ,i a s (M z ) . 



(6) 



The parameters c„ , are fitted to the theoretical cross sec- 
tions obtained for different choices of a s (Mz), where the 
range a s (Mz) — 0.1 15 to 0. 123 is being used. The cross 
section measured in bin i is mapped to the parametrisa- 
tion such that the obtained value of a s (M z ) can be ob- 
tained easily. The experimental uncertainty on a s (Mz) 
is obtained by projecting the uncertainty of the mea- 
sured cross section value on the parametrisation. A 
combined value of a s (Mz) from N bins is obtained by a 
^ 2 -fit. The correlated systematic uncertainties are taken 
into account by re-fitting a s (Mz) for all bins with the 
systematic shifts applied to the measured cross section 
values [17|. 

The theoretical uncertainties due to missing higher 
orders are estimated by the band method suggested by 
Jones et al. IPT81 . In this method the theoretical predic- 
tions are recalculated using the obtained central value 
of a s (Mz) with the renormalisation and factorisation 
scales varied up and down by a factor of two. The ob- 
tained differences in the prediction are compared to the 
nominal theory with varying a s (M z ). The correspond- 
ing smallest and largest values of a s (Mz) are then stored 
for each bin. In a last step the validity of the resulting 
values of a s (Mz) is tested for each measurement bin and 
only the allowed values of a s (Mz) are used to set the 
theoretical uncertainty on a s {M z ). The uncertainty on 
a s (Mz) due to the parametrisation of the proton PDFs is 
obtained in the same way. 

5.3. Comparison of both methods 

The same NLO calculations are used for the a s deter- 
minations performed by both experiments. For the DIS 
case calculations based on the Catani-Seymour dipole 
subtraction [19] are employed, as implemented in the 
program NLOJet++ J20] |2TJ. In the case of photo- 
production, calculations using the phase-space-slicing 
method lt22l as implemented in the program by Klasen, 
Kleinwort and Kramer l23l are used. 
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However, as outlined above, the exact methods used 
by the two experiments for the determination of a s (Mz) 
are quite different. The HI -type fit takes the systematic 
uncertainties into account during the fit through the in- 
troduction of penalty terms, which allows a cross-check 
of the determined uncertainties due to systematic mis- 
measurements. This is not done in the ZEUS method. 
On the other hand, in the ZEUS case the correlation be- 
tween a s and the gluon density is correctly taken into 
account, which is neglected in the case of HI. In the 
latter case the assumption on the gluon PDF can only 
be tested after the fit by using PDF sets obtained with 
different values of a s (Mz). 

The theoretical uncertainties are determined in the 
ZEUS fit based on theory predictions only, whereas in 
the offset method employed by HI the experimental un- 
certainties enter the fit when estimating them. Together 
with the band method this results in a smaller theoretical 
uncertainty by about 30-50% 1 16 1 in the ZEUS case. 

6. Experimental environment 

6.1. The HERA collider 

The HERA collider was located in Hamburg, Ger- 
many and has been in operation between 1992 and 2007. 
After 15 years of successful operation the two colliding 
beam experiments HI and ZEUS each have collected 
about 0.5 fbr 1 of ep scattering data. Between the years 
2001-2003 a machine upgrade has been carried out, 
which brought an increase of the instantaneous luminos- 
ity by a factor of 4—5 as well as the possibility of longi- 
tudinally polarised beams. This shutdown also has been 
used to replace and upgrade some of the sub-detectors 
of the experiments. At the end of the HERA-1 phase 
as well as during the HERA-2 phase the beam energies 
were 27.6 GeV for the electrorQbeam and 920 GeV for 
the proton beam, resulting in a centre-of-mass energy of 
about 318 GeV. 

6.2. The HI and ZEUS experiments 

The two colliding-beam experiments HI and ZEUS 
were multi-purpose particle detectors with an asymmet- 
ric instrumentation, which reflected the fact that the 
hadronic centre-of mass system was boosted in the pro- 
ton direction. Both were equipped with silicon strip de- 
tectors close to the interaction point and large jet drift 



'The HERA collider was operated with electrons and positrons, 
but since the charge of the lepton beam is of no importance for the 
physics discussed here, the term electron will be used generically for 
electrons and positrons. 



chambers inside a magnetic field of 1.16 and 1.43 T, for 
the HI and ZEUS case, respectively. The magnetic field 
was provided by superconducting coils in both experi- 
ments, with the difference being the placement of them. 
In the case of HI the superconducting coil surrounded 
the calorimeter stacks, while in the ZEUS case it was 
placed between the tracking detectors and the calorime- 
ters. However, the largest difference between the two 
detectors were the calorimeters employed. HI used a 
high-granularity liquid-argon calorimeter with lead and 
steel plates as absorbers, whereas ZEUS was equipped 
with a compensating calorimeter based on depleted ura- 
nium as absorber and scintillating fibres for the signal 
detection. This difference is reflected in the achieved 
energy resolutions. While HI achieved a better resolu- 
tion for electrons with 11%/Vx? compared to 18%/Vc 
in the ZEUS case, for hadrons ZEUS achieved a bet- 
ter resolution with 35%/ Ve compared to 55%/ Ve for 
HI, with all resolutions obtained from test beam data. In 
both experiments large area muon chambers surrounded 
the calorimeters for the rejection of cosmic ray back- 
ground, muon identification and the measurement of en- 
ergy leakage from the calorimeters. 

7. Event reconstruction and jet measurement 

In DIS the kinematics can be completely determined 
by the measurement of the scattered electron. There- 
fore, the electron identification and reconstruction is 
crucial for NC measurements. The HI and ZEUS col- 
laborations have achieved excellent precision of the 
electron energy measurement with experimental uncer- 
tainties of the order of 1%. Also for jet measurements 
in DIS the electron reconstruction is crucial, since the 
scattered electron is used for the reconstruction of the 
event kinematics and the determination of the boost to 
the Breit frame. 

In jet measurements the dominating experimental un- 
certainty is usually the jet energy scale uncertainty due 
to the steeply falling cross section as function of jet Pj. 
In the following the jet reconstruction for the two exper- 
iments is reviewed. 

7.1. Jet reconstruction at ZEUS 

The compensating calorimeter of the ZEUS experi- 
ment with its high resolution for hadronic showers was 
ideally suited for jet measurements. However, due to 
the material from the beam pipe, tracking system and 
superconducting coil in front of the calorimeter on aver- 
age about 20% of the jet's transverse energy was lost 
and needed to be corrected for Il24l . This has been 
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Figure 2: Ratio of the transverse momentum of the calibrated 
hadronic final state Pj to the reference measurement P^ 1 as function 
of the jet energy (a) and the pseudorapidity (b) as obtained by the HI 
collaboration. The ratio of the data to MC comparison is shown at the 
bottom of each plot. 
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Figure 3: The measured cross section for inclusive dijet (a), (b) and 
trijet production (c), (d) at low Q 2 as function of Q 2 and the average 
transverse momentum (Pj) by the HI collaboration. The measure- 
ments are compared to NLO calculations corrected for hadronisation 
effects. 



achieved by deriving dead-material corrections from 
Monte Carlo (MC) simulations in photoproduction sam- 
ples using calorimetric and tracking information |25|. 
Residual corrections were obtained by balancing the 
scattered electron with the hadronic final state in NC 
events, individually for data and simulated event sam- 
ples l26l . With these corrections applied, a jet energy 
scale uncertainty of 1% for jets with Pj > 10 GeV has 
been achieved for ZEUS already in the year 2002. The 
recent jet measurements by the ZEUS collaboration de- 
scribed in this document use only calorimetric informa- 
tion for the jet reconstruction. 

7.2. Jet reconstruction at HI 

Due to the non-compensating nature of the HI 
calorimeter the energy calibration is more involved and 
only recently a jet energy scale uncertainty of 1 % has 
been achieved. The energy calibration procedure is per- 
formed in several steps. 

In the first step a software weighting is applied, which 
corrects for the large fluctuations of the invisible en- 
ergy fraction in hadronic showers. This ;r -weighting 
improves the energy reconstruction and resolution by 
weighting energy deposits depending on the energy den- 
sity in single calorimeter cells |[27|. However, the cali- 
bration coefficients for this method were obtained from 
MC simulations where the absolute energy scale could 
not be sufficiently well simulated. This resulted in dif- 
ferences between data and MC simulations of about 4%. 

In order to improve the energy reconstruction an 
energy -flow algorithm was developed [28 , 29 1. This al- 
gorithm avoids double-counting of energy by compar- 
ing the momentum resolution of a reconstructed track 



with the expected resolution of a matching energy de- 
posit in the calorimeter. If the resolution of the recon- 
structed track is better, an amount of energy compati- 
ble with the track measurement is subtracted from the 
calorimetric measurement. Otherwise the track is dis- 
carded and the calorimetric measurement is kept. 

Recently this algorithm has been improved by sep- 
arating calorimetric energy deposits originating from 
electromagnetic and hadronic showers with the help of 
neural networks 113 PL Based on the obtained electro- 
magnetic probability a calibration method has been de- 
veloped which can be applied independently to data and 
simulated events. Calibration constants for individual 
calorimeter clusters have been obtained as a function 
of their energy and angular region 1 30 1 . A comparison 
of the obtained transverse momentum of the calibrated 
hadronic final state with a reference measurement is 
shown in Fig.[2]for data and simulated events PP . An 
absolute energy reconstruction of 2% and an agreement 
between data and MC simulation of 1% is achieved over 
the full range of accessible jet energies and pseudora- 
pidities. 

8. Jet production at low Q 2 

In a recent analysis by the HI collaboration jet pro- 
duction in the Breit frame is measured at low virtuality 
of the exchanged boson, 5 < Q 2 < 100 GeV 2 Il32l . 
In this analysis data from the HERA-1 running phase 
is used, however, at low Q 2 the cross section is large 
such that even at high jet Pj the statistical uncertain- 
ties are small. The high statistics allowed for a double- 
differential measurement of inclusive jet, dijet and tri- 
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jet production as function of Q 2 and jet Pj (the av- 
erage transverse momentum (Pj) in case of the dijet 
and trijet measurements). An invariant mass cut of 
Mi 2 > 18 GeV for the two jets with highest P T in the 
event is required for the dijet and trijet measurements 
in order to obtain reliable predictions from NLO calcu- 
lations. The measured dijet and trijet cross sections as 
function Q 2 and (Pj) in the Breit frame are compared to 
NLO calculations in Figure [3] The turnover of the tri- 
jet distribution at (P T ) around 12 GeV is due to the cut 
on the invariant mass and is well described by the NLO 
calculations. The NLO calculations are obtained with 
the choice yu r = -<J((Pt} 2 + Q 2 )/2. As an alternative, the 
choice fi r = (Pj) has been tested but it is disfavoured by 
the data in regions where Q is larger than (Pj). The the- 
oretical uncertainties, mostly due to the variation of fi T , 
are large in the phase space of this measurement. They 
are of the order of 30% at low (Pj) and about 10% at 
high (Pi), whereas the total experimental uncertainties 
are about 6-10%. 

As a first step, values of a s (Mz) are extracted for each 
of the 62 measured cross section points for inclusive jet, 
dijet and trijet production individually. Good agreement 
between the obtained values is observed. In a second 
step a combined fit to all 62 data points is performed 
using Eq. Q, resulting in a value of 

a s (M z ) = 0.1160 ± 0.0014 (exp.) 

+ 0.0016 (pdf) (theo.), 

with a good fit quality of^- 2 /ndf = 49.8/61. The large 
theoretical uncertainty is a consequence of the large un- 
certainty due to missing higher orders in the NLO calcu- 
lation. Smaller theoretical uncertainties can be obtained 
by studying jet production at high Q 2 . 

9. Jet production at high Q 2 

The HI and ZEUS collaborations have measured jet 
production for Q 2 values larger than 150 and 125 GeV 2 , 
respectively. Similarly as in the low Q 2 case, inclu- 
sive jet, dijet and trijet production cross sections have 
been measured with high precision, made possible by 
the larger data samples available from the HERA-2 run- 
ning phase. 

9.1. Inclusive jet production 

Preliminary results on inclusive jet production at high 
Q 2 using HERA-2 data have been released by the HI 
and ZEUS collaborations ll3Tl[33ll . While the range cov- 
ered in Q 2 is similar in both analyses, the additional re- 
quirements on the phase space in which the measure- 
ments are performed are quite different. While the HI 
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Figure 4: The ratio of NLO calculations using different proton PDFs 
to the measured cross sections of inclusive jet production at high Q 2 
from HI. The theoretical uncertainty due to missing higher orders, 
estimated by a variation of p. t and pf, is shown for the calculation 
using HERAPDF 1.5 only. 



collaboration restricts the inelasticity to 0.2 < y < 0.7 
and the pseudorapidity range of the jets as measured 
in the laboratory rest frame to -1.0 < rji^ < 2.5, in 
the ZEUS analysis the measurement is performed for 
| cosy;, | < 0.65 and the jet pseudorapidity requirement 
in the Breit frame -2.0 < ?7Breit < L5. The variable jh is 
the scattering angle of the hadronic final state and cor- 
responds to the polar angle of the scattered quark in the 
LO picture of DIS. These requirements are imposed in 
order to have an experimentally well controlled boost 
to the Breit frame and jets well within the geometri- 
cal detector acceptance. Both collaborations measured 
the cross section for inclusive jet production double- 
differentially as function of Q 2 and jet Pj. 

The ratios of the calculated cross sections at NLO 
QCD, obtained with different proton PDFs, to the mea- 
sured ones by HI are shown in figure|4] The NLO calcu- 
lations give a good description of the data. The largest 
spread of the predictions obtained with different PDF 
sets is observed for Pj > 30 GeV for the Q 2 bins 4 and 
5, corresponding to 400 < Q 2 < 5000 GeV 2 . In this 
region inclusive jet production is dominated by the va- 
lence quark distributions at 0.1 < x p < 0.7, where x p is 
the longitudinal momentum fraction of the proton taken 
by the valence quark, as used by the NLO calculations. 
Since at small values of Q 2 and Pj the cross section is 
dominated by gluon-induced processes, the data have 
the potential to constrain the valence quark at high x p 
and provide direct sensitivity to the gluon distribution 
and a s . 

The extraction of a s (M z ) by HI is performed for all 
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24 measured bins, and the resulting value is 



ZEUS 



a s (M z ) = 0.1 190 + 0.0021 (exp.) 

± 0.0020 (pdf) 



+0.0050 
-0.0056 



(theo.) . 



where the theoretical uncertainty is smaller by nearly 
a factor of two with respect to the low Q 2 measure- 
ment. Further increasing the Q 2 cut for the determi- 
nation of a s (M z ) improves the theory uncertainty with 
the drawback of a larger experimental uncertainty. In 
order to minimise the total uncertainty on a s , the ZEUS 
collaboration performed a determination of a s (Mz) for 
Q 2 > 500 GeV 2 . The fit is made using four data points 
of the single-differential do-j el /dQ 2 measurement. The 
resulting value of a s (M z ) is 



a s (M z ) = 0.1208 !^ (exp.) 

± 0.0008 (pdf) 



: 0.0022 (theo.) , 



where the small uncertainty due to the parametrisation 
of the proton PDFs is obtained using the ZEUS-S 04l 
PDF set. The theoretical uncertainty, mostly due to 
missing higher orders but also including uncertainties 
of the hadronisation corrections, is only half as large as 
the corresponding value obtained by HI. 

9.2. Dijet production 

Although the cross section of inclusive dijet produc- 
tion in the Breit frame is calculated using the same di- 
agrams as inclusive jet production, there are important 
differences. While in inclusive jet production all jets in 
a given pseudorapidity region with Pj above a certain 
threshold contribute, an event contributes to the dijet 
cross section only if at least two jets above a minimum 
Pt are found within the given pseudorapidity range. 
This has important consequences. In the inclusive jet 
case events contribute to the cross section where only 
one jet lies in the central region and is balanced in trans- 
verse momentum by one or more jets with large pseudo- 
rapidity. These configurations have larger higher-order 
corrections, which is reflected in a larger theoretical un- 
certainty due to the variation of yU r for inclusive jets. In 
the dijet case either an invariant mass or an asymmetric 
P T cut for the two leading jets^Jis required if the cross 
section obtained from a fixed-order calculation should 
be reliable. In the case of a symmetric Pt cut with- 
out a restriction on the invariant mass M\i, there are 
regions in phase space sensitive to soft gluon emission 



The jet with the highest Pt in the event is referred to as the lead- 
ing jet. 
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Figure 5: The relative differences between the measured double dif- 
ferential dijet cross sections by ZEUS and NLO calculations with dif- 
ferent choices of p. Y . The theoretical uncertainty, shown as dashed 
area, is only shown for the central prediction with p Y = y]Q- + (Pt) 2 - 



and all-order resummations would be needed [ 13 1. With 
either an asymmetric Pt cut or a requirement on M\i 
fixed-order calculations have been shown to describe 
dijet data sufficiently. However, a smaller uncertainty 
due to missing higher orders can be obtained with a cut 
on the invariant mass. Therefore, the most recent di- 
jet measurements from HI and ZEUS define the dijet 
phase space through the requirements Pt > 5 GeV and 
M12 > 16 GeV in the case of the HI analysis BTl and 
P T > 8 GeV and M 12 > 20 GeV in the case of the ZEUS 
analysis 11351 . 

The relative differences between dijet data measured 
double differentially by the ZEUS collaboration and 
NLO calculations with different choices of the scale 
Hr are shown in figure Fig. [5] A scale choice of 
A<r = VS 2 + (^t) 2 , identical to the choice by HI up to 
a factor of 1/V2, leads to a good description of the 
data over the full region of Q 2 and (Pt). The choice 
of ju r = (Pt) fails to describe the data for large values of 
Q 2 when at the same time (P T ) is small. 

As mentioned above, the theoretical uncertainties due 
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Figure 6: Comparison of the trijet cross sections measured by the HI 
collaboration with NLO calculations using HERAPDF 1.5, The cross 
sections are shown as function of (Pj) in different Q 2 bins. The ratio 
of theoretical predictions to data is shown at the bottom of each plot. 



to missing higher orders are smaller in the dijet analysis 
than in the inclusive jet case. The HI collaboration de- 
termined a s (M z ) from a dijet measurement with a si- 
multaneous fit to 24 bins and obtained a value of 

a s (M z ) = 0.1 146 ± 0.0022 (exp.) 

± 0.0021 (pdf) !°;°°£(thea.), 

where the theoretical uncertainty is about 20% smaller 
than the one obtained for the inclusive jet measurement. 

9.3. Trijet production 

The improved reconstruction of the hadronic final 
state and jet energy scale, together with the high statis- 
tics from the HERA-2 running, allowed for the first time 
a double-differential measurement of trijet production at 
high Q 2 1311 . In this analysis the same requirement on 
the invariant mass of the two leading jets as in the dijet 
analysis is required, Mn > 16 GeV, and all jets are re- 
quired to have P T > 5 GeV. Because of this, the events 



used to measure the trijet cross sections are a subsample 
of the events used for the dijet analysis. This is benefi- 
cial for future QCD analyses where both measurements 
are used. 

The measured trijet cross sections have larger ex- 
perimental uncertainties than the inclusive jet or dijet 
cross sections, due to larger model uncertainties and a 
larger effect of the jet energy scale uncertainty. The 
measurement has total experimental uncertainties be- 
tween between 7% at low values of (Pj) and Q 2 and 
15% at high values of (Pj) and Q 2 . Also the theoret- 
ical uncertainties are larger than in the dijet case by 
about 50% and in most bins twice as large as the ex- 
perimental uncertainties. The NLO calculations provide 
a good description of the trijet data with the choice of 
p x = ~<J((Pt} 2 + Q 2 )/2. A comparison of the measured 
trijet cross sections with NLO calculations using a value 
of a s (M z ) = 0.118 is shown in Fig. [6] 

The obtained value of a s (M z ) from trijet cross sec- 
tions has the smallest experimental and PDF uncertainty 
when compared to the values obtained from the inclu- 
sive jet and dijet measurements by HI. This is due to 
the fact that the trijet cross section is already in LO pro- 
portional to c". The value of a s (M z ) obtained is 

a s (M z ) = 0.1196 ± 0.0016 (exp.) 

±0.0010 (pdf) !°;°°|(theo.), 

being in good agreement with the values obtained from 
the inclusive and dijet analyses. 

9.4. Normalised jet cross sections 

One way to improve the experimental precision fur- 
ther is to measure ratios of cross sections. For jet pro- 
duction in NC DIS the ratio of jet cross sections to the 
inclusive NC cross section is an obvious choice. These 
normalised jet cross sections can be understood as jet 
rates in DIS. The advantage of these observables lies 
in the cancellation of some experimental uncertainties 
such as the trigger uncertainty and the uncertainty of 
the luminosity measurement. Some other uncertainties 
such as the model dependence of the acceptance correc- 
tion or the uncertainties related to the electron measure- 
ment may cancel partially in the ratio, leading to a fur- 
ther improvement of the experimental precision. This 
has been impressively demonstrated by a measurement 
of normalised inclusive jet, dijet and trijet cross sec- 
tions [16 1 by the HI collaboration. A simultaneous fit 
to all 54 data points gives 

a s (M z ) = 0.1 168 ± 0.0007 (exp.) 

±0.0016 (pdf) ™6(theo.), 
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with a good^- 2 /ndf of 65/53. The achieved experimen- 
tal precision of 0.6% raises the hope of an uncertainty of 
0(1%) once next-to-next-to-leading order (NNLO) cal- 
culations become available. 

10. Jets in photoproduction 

With the study of jet production in y*p interactions 
it is possible to gain information on the partonic con- 
tent of the photon in addition to a s and the structure 
of the proton. In contrast to DIS, there is only one 
hard scale involved since Q 2 ~ GeV 2 in photopro- 
duction. The only available choice for the scales yu r 
and //f in the NLO calculations is therefore the jet Pj. 
Although this simplifies the calculations at first sight, 
the hadronic structure of the photon makes y*p interac- 
tions similar to hadron-hadron interactions which com- 
plicates matters. This also means that additional per- 
turbative and non-perturbative effects such as multiple- 
interactions can play a role in the calculations. How- 
ever, the advantage of jet production in y*p compared 
to hadon-hadron interactions is the absence of multiple- 
interactions for the point-like contribution. This allows 
the study of multiple-interaction models in an interest- 
ing transition region. 

Recently the ZEUS collaboration measured inclusive 
jet production in photoproduction [36 1, using the full 
HERA-2 data. In this analysis the selected events are 
restricted to Q 2 < 1 GeV 2 , and the yp centre-of-mass 
energies are in the range of 142 < W yp < 293 GeV, 
which corresponds to 0.2 < y < 0.85. The selected jets 
are required to have Pj > 17 GeV and have to lie in 
the pseudorapidity region -1.0 < 77^ < 2.5. The NLO 
calculations are corrected for hadronisation effects by 
using the ratio of cross sections on the level of stable 
hadrons to the partonic cross sections using the Pythia 
and Herwig MC. Disagreements of up to 40% between 
the data and the NLO predictions can be found at large 
values of pseudorapidity, rj > 2, and small values of jet 
transverse momenta, Pj < 30 GeV l36l . This region 
corresponds to small values of jt° bs , see Eq. (|3j, where 
the resolved photon is expected to give a large con- 
tribution to the cross section and in addition multiple- 
interactions may occur. In order to study the differences 
between the data and the predictions in the context of 
multiple-interactions, an additional correction factor is 
calculated. It is computed as the ratio of cross sections 
on the level of stable hadrons using Pythia samples with 
multiple-interactions and samples without them. The 
probability of a secondary interaction to occur depends 
strongly on its allowed value for the minimum trans- 
verse momentum n!? c . . Three sets of correction fac- 



ZEUS 




Figure 7: The inclusive jet cross sections in photoproduction as 
function of jet pseudorapidity ^ ct . The measured cross sections are 
compared to NLO calculations corrected for hadronisation effects and 
for non-perturbative effects from multiple-interactions (NP). The lat- 
ter corrections were obtained with Pythia using different assumptions 
on the minimum Pj of the secondary interaction. 

tors are obtained for p55 c . = 1, 1.5 and 2 GeV. The 

1 i,min ' 

comparison of the modified NLO calculations with the 
data is shown in Fig. [7] While the NLO calculations 
without additional corrections lie significantly below 
the data, with the disagreement increasing as function 
of rj, the calculations with an additional correction using 
reL ec =1.5 show excellent agreement with the mea- 
sured cross sections. However, also a different choice 
of photon PDFs can reduce the discrepancy between 
the data and the NLO calculations. Increasing the re- 
quirement on the minimum jet Pj reduces the influence 
of hadronic component of the photon and therefore re- 
duces the contribution of multiple-interactions and the 
uncertainty due to the photon PDFs. This behaviour 
is confirmed in this analysis, where a requirement of 
Pj > 21 GeV for jets results in reasonable agreement 
between the NLO calculations corrected for hadronisa- 
tion effects only and the data, as shown in Fig. [8] The 
total theoretical uncertainty shown as hatched area con- 
sists of the uncertainties due to the scale variations, the 
proton and photon PDFs, and the hadronisation correc- 
tions, added in quadrature. It is larger than the total 
theoretical uncertainty in the case of the high Q 2 mea- 
surements due to the uncertainty from the photon PDFs, 
which can be as large as 13% at high values of rj. 
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The value of a s (Mz) is extracted then from the 
dcr/dP T distribution for 21 < P T < 71 GeV, where the 
upper cut on Pj is motivated by a relatively large uncer- 
tainty due to the proton PDFs for higher Pj values. The 
value obtained is 



using the GRV-HO 



set of photon PDFs. At 



a s (M z ) = 0.1206 



+0.0023 
-0.0022 



(exp.) 



: 0.0030 (pdf) ^(th.h 



where the uncertainties due to the proton and photon 
PDFs have been added in quadrature. The uncertainties 
of the photon PDFs have a much larger effect on a s (Mz) 
than the uncertainties of the proton PDFs. The rela- 
tive uncertainty on a s (Mz) related to the photon PDFs 
is 2.3%, which has to be compared with 1% due to the 
proton PDFs. 

A possibility to further constrain the photon PDFs 
and reduce their uncertainties is dijet production in pho- 
toproduction, where it is possible to calculate jc° bs as de- 
fined in Eq. (|3). In a new analysis using HERA-2 data, 
the ZEUS collaboration has measured dijets in y*p in- 
teractions 1 37] with an asymmetric P T requirement of 
P T ,i > 21 GeV and P T ,2 > 17 GeV. The measured dijet 
cross sections as function of x° bs are compared to the 
NLO calculation in Fig. [9] The highest jet Pj in the 
event is chosen for the scales p r and fif in the NLO cal- 
culations and a good description of the data is observed 



small values of x° bs , as the contribution of the resolved 
photon becomes more important, a large spread of the 
predictions using different photon PDF sets is observed. 
Dijet data from HERA may therefore be important for 
future photon PDF determinations and can be used for 
direct tests of the gluonic content of the photon. 



11. Summary of recent a s values from HERA 

A stringent test of QCD can be performed by testing 
the running of the coupling a s at different scales yU, as 
given by Eq. ([T]). Jet data from HERA provide a good 
possibility for this test, since they cover a large range of 
scales in a single process. Tests of the running of a s are 
obtained by directly fitting a s (fi T ), instead of evolving 
a s to the value of M z . Results from such a fit are shown 
in Fig. 



10 



where a s (p Y ) has been fitted to normalised 
jet data at high Q 2 ifTBI and absolute jet cross sections 
at low Q 2 [ 32 1 . For this fit the double differential jet 
cross sections d 2 cr/dPjdQ 2 have been used, where the 
values of a s obtained in the different Q 2 bins have been 
combined. Very good agreement between the obtained 
values of a s (jj. r ) with the predicted running of a s can be 
observed, where the data cover scales from about 6 to 
65 GeV. Interestingly, the a s values obtained in the low 
Q 2 analysis agree with the extrapolated a s (ju r ) from the 
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Figure 10: Running of a s (j*r) obtained from multi-jet cross sections 
in DIS ( 32 ] . The theoretical uncertainty, mostly due to missing higher 
orders, has been extrapolated using the two-loop solution of the renor- 
malisation group equation from the high Q 2 analysis to the small Q 2 
analysis. This explains the relatively small theoretical uncertainty at 
small values of p c . 



high Q 2 analysis within the experimental uncertainties. 
This is quite striking since the theoretical uncertainties 
in the low Q 2 analysis are twice as large as the ones 
in the high Q 2 analysis, where only the latter ones are 



shown as grey band in Fig. 10 



A similar study ll37l has been carried out by the 
ZEUS collaboration, where the running of a s (ju r ) is 
tested with inclusive jet data in photoproduction and 
high Q 2 NC DIS. The result is shown in Fig [TT] In 
this case, the fit has been performed using the measured 
dcr/dPj distributions, and the average jet Pj measured 
in each bin has been chosen for the scale jj r . Good 
agreement between the obtained values of a s (fj. r ) and the 
two-loop solution of the renormalisation group equation 
is observed also in this analysis. The transition between 
photoproduction and DIS is smooth and does not show 
any peculiarities. 

A summary of all values of a s (M z ) given in this 
document is shown in Fig. 12 The obtained values 



show good agreement with the world average [40]. 
The HERA measurements of a s (Mz) from jet data 
are clearly dominated by the theoretical uncertain- 
ties, which are mostly due to missing higher orders. 
However, the experimental uncertainties are compara- 
ble with the precision of the world average and once 
NNLO calculations become available, determinations 
of a s (Mz) with a precision at the percent level will be 
feasible using jet data from HERA. It should be noted 
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Figure 11: Running of a s (ji r ) obtained from inclusive jet cross sec- 
tions in photoproduction and NC DIS f37l. The experimental uncer- 
tainties are shown with their uncorrelated and correlated parts. The 
solid line corresponds to the two-loop solution of the renormalisation 
group equation with a s (Mz) = 0.1206. 



that the precision of the world average is mostly driven 
by lattice calculations and the analysis of r-decays. 
When comparing the precision of the recent a s deter- 
minations from HERA jet data with other collider data, 
the achieved precisions are competitive. 

12. Conclusion 

The last few years have brought a large improvement 
in terms of experimental precision for jet measurements 
from HERA. This is partly due to the higher statistics of 
the HERA-2 data, but mostly because of achievements 
of the experimental collaborations in terms of their re- 
construction and calibration algorithms. Both, the HI 
and ZEUS collaborations have reached the goal of a 
1% jet energy scale uncertainty. They are now in the 
process of finalising their physics programmes which 
include high precision jet measurements, with the first 
ones published recently. 

These measurements include analyses of inclusive 
jet, dijet and trijet production in different kinematic 
regimes and provide stringent tests of pQCD. Even 
more importantly, these data will be important for future 
QCD analyses. The potential of these data can already 
be seen by the precision determinations of a s (Mz) with 
experimental uncertainties of 0(1%). However, the the- 
oretical uncertainties, which are mostly due to missing 
higher orders, are limiting the precision of the extracted 
values of a s (Mz) to about 3-4%. The obtained values 
of a„(Mz) from the different analyses show excellent 
agreement. This is quite remarkable considering that 
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Figure 12: Comparison of the most recent values of a s (Mz) obtained 
from multi-jet production cross sections at HERA and the world aver- 
age (40). 

they have been obtained using different methods and 
data in very different kinematic regimes from two dif- 
ferent experiments. 

This opens the possibility for a determination of 
®s(Mz) using the latest jet data from HI and ZEUS 
simultaneously, which may improve the experimental 
precision further. The last attempt in this direction has 
been made in 2007 and showed a reduction of the ex- 
perimental uncertainty of about 30% with respect to the 
individual determinations [41 1. 

Another interesting possibility presents itself by a 
combination of jet cross sections in a similar way as 
the combination of inclusive DIS data [4|. There the 
final cross sections are obtained by a fit to the individ- 
ual data sets. In this way not only the statistical uncer- 
tainties improve, but also the systematic uncertainties 
are potentially reduced because the two experiments are 
allowed to cross-calibrate one another. This approach 
would have the advantage of providing jet data with the 
ultimate precision, which can be used in future QCD 
analyses with NNLO accuracy. 
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